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ABSTRACT 

We investigate the impact of an early population of massive stars on their 
surroundings. Dissociation of molecular hydrogen by strong UV emission from 
such stars is expected to produce a global transition in the cooling mechanism 
of minihalos at a redshift of approximately 30, strongly inhibiting star formation 
until more massive halos can collapse. Furthermore, chemical enrichment from 
Pop III supernovae produces a second transition at z ~ 15 — 20, when the mean 
metallicity of the universe exceeds a critical threshold and Pop III star formation 
gives way to Pop II. We show that observations of high redshift supernovae with 
the Next Generation Space Telescope (NGST) have the potential to trace the 
cosmic star formation rate out to z > 20, provided that Pop III supernovae 
are at least as bright as, and ideally brighter than, type la supernovae. We 
also propose a mechanism for the formation of a novel population of extremely 
low metallicity stars of intermediate mass at very high redshifts, which we term 
Pop II. 5. In our model shock compression, heating, and subsequent cooling to 
high density reduces the fragment mass in primordial gas to ~ 10 M Q , allowing 
low mass stars to form. We predict the number density of relic Pop II. 5 stars in 
the Milky Way halo today and find that, with certain assumptions, there should 
be ~ 10 kpc~ 3 in the solar neighborhood. 

Subject headings: cosmology: theory — early universe — galaxies: formation - 
intergalactic medium — stars: formation 

1. Introduction 

An important challenge in modern cosmology is to understand how the cosmic "dark 
ages" ended (for a recent review, see Barkana & Loeb 2001). There is growing theoretical 
evidence indicating that the first luminous objects to form in the universe were very massive 
stars with typical masses M* > 100 M (Bromm, Coppi & Larson 1999, 2002; Abel, Bryan & 
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Norman 2000, 2002; Nakamura & Umemura 2001). These stars formed out of metal-free gas 
in dark matter (DM) halos of mass ~ 10 5 - 10 6 M at redshifts z > 20 (e.g., Tegmark et al. 
1997; Fuller & Couchman 2000). Typically, simulations predict that the halos hosting this 
process should contain either one or a few dense massive clumps of baryonic matter, with the 
most massive clump situated nearest to the center of the halo being the one to collapse first 
and to presumably form a very massive star (VMS). On the other hand, conditions in the 
universe less than 10 9 yr after the big bang must have allowed for the formation of low-mass 
stars, with typical masses of ~ 1 M , as is implied by the ages of the oldest globular clusters 
in the Galactic halo with metallicities of Z ~ 10~ 2 Z Q (e.g., Ashman & Zepf 1998; Bromm & 
Clarke 2002). In this paper, we investigate the question: How and when did the transition 
from the formation of very massive stars at high redshifts to that of more normal, low-mass 
stars at later times take place? 

Observations have recently provided hints to the character of star formation at z > 5. 
The abundance of C IV in the low-column density Lya forest indicates that the intergalactic 
medium (IGM) was enriched with heavy elements to a level of Z ~ 1O _3 ' 5 Z already at 
z ~ 5 (Songaila 2001). In addition, abundance patterns in Lja systems out to z ~ 4.6 are 
suggestive of an early, prompt nucleosynthetic inventory of heavy elements from a generation 
of very massive stars (Qian, Sargent & Wasserburg 2002). 

Elucidating the cosmic star formation history in the high redshift universe is crucial 
for predicting what observatories such as the Next Generation Space Telescope (NGST) will 
observe at z > 5 less than a decade from now. In particular, it is important to ascertain the 
rates and properties of high redshift supernovae (SNe) which are likely to be the brightest 
beacons heralding the end of the "dark ages" (Miralda-Escude & Rees 1997). Similarly, 
gamma-ray bursts (GRBs) are expected to trace the formation history of massive stars out 
to very high redshifts (Totani 1997; Wijers et al. 1998; Blain & Natarajan 2000). In fact, 
the top-heavy initial mass function (IMF) predicted for the first stars favors the massive 
stars which are the likely source of GRB progenitors. Constraining how stars form at high 
z, therefore, is of great relevance for interpreting the results of the upcoming Swift satellite 
(Bromm & Loeb 2002). 

Here, we examine processes related to star formation in the high redshift universe using 
a simple analytical model that synthesizes a number of recent theoretical and observational 
results into a coherent framework. This idealized model is intended to emphasize the basic 
physical mechanisms, and we will complement it with detailed numerical simulations in 
future work. 

The organization of the paper is as follows. In §2 we describe our high-redshift star 
formation model. A discussion of the various stellar populations that form at high z is given 
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in §3, while §4 presents observational consequences. Finally, §5 contains our conclusions and 
the outlook for future work. 



In general, VMSs will impart three forms of feedback on the IGM (e.g., Ciardi et al. 
2000b): radiative, chemical, and mechanical. Owing to their high effective temperatures, 
T e fj rsj 10 5 K, these stars produce a large quantity of UV photons, roughly 10 times as many 
per stellar baryon as those from a "normal" stellar population characterized by a Salpeter 
IMF (Tumlinson k Shull 2000; Bromm, Kudritzki & Loeb 2001b). Upon its death after only 
~ 3 x 10 6 yr, a VMS explodes as a supernova and enriches the IGM with heavy elements, 
provided the star is not more massive than ~ 260 M Q , in which case it is predicted to 
collapse entirely into a massive black hole without any concomitant metal ejection (Fryer, 
Woosley & Heger 2001). We propose that this radiative and chemical feedback gives rise to 
two fundamental transitions in the cosmic star formation history at high z. The feedback 
from mechanical energy input into the IGM might lead to a distinct new population of stars, 
to be discussed in §3. 

To assess the global feedback exerted by VMSs on the IGM, we need to know their 
formation rate. We calculate this rate according to the "single-burst" model of VMS for- 
mation introduced by Santos, Bromm & Kamionkowski (2002). This model adopts the 
observationally favored scenario in which the formation of cosmic structure has progressed 
hierarchically from small to large scales as described by the cold dark matter (CDM) model 
with a cosmological constant. Specifically, we assume a ACDM model with density param- 
eters in matter Q m = 1 — Q\ = 0.3 and in baryons Qb = 0.045, a Hubble parameter of 
h = if /(100kms _1 Mpc -1 ) = 0.7, and a scale-invariant power spectrum of density fluctua- 
tions with a normalization og = 0.9. The abundance and merger history of the DM halos 
is given by the extended Press-Schechter formalism (Lacey & Cole 1993). Within merging 
halos, stars form in gas that is able to cool efficiently, via either molecular or atomic hydro- 
gen lines. These cooling mechanisms become effective at different minimum temperatures, 
^crit ~ 400 K and ~ 10 4 K for molecular and atomic hydrogen, respectively (e.g., Barkana & 
Loeb 2001). These temperature thresholds correspond to redshift- dependent minimum DM 
halo masses 



where \i is the mean molecular weight (/i = 0.6 for ionized gas, and \i = 1.2 for neutral gas). 



2. Epochs of Star Formation at High Redshift 



2.1. Cosmic Star Formation History 
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Fig. 1. — History of comoving star formation rate (SFR) density in units of M yr _1 Mpc _1 
as a function of redshift. Solid line: Formation rate of VMSs in halos that cool due to atomic 
hydrogen. Dotted line: VMS formation rate in halos that cool due to molecular hydrogen. 
In calculating these rates, VMS formation was assumed to proceed in a single-burst mode 
with an efficiency of i]* = 10% (see Santos et al. 2002). Dashed line: Analytic fit to the 
star formation history derived from cosmological simulations that incorporate a multi-phase 
model of quiescent star formation (from Springel & Hernquist 2002). 
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Finally, the star formation efficiency is chosen to be 77* = O.fO. The "single-burst" model 
now assumes that VMS formation is accompanied by strong negative feedback effects on the 
star- forming gas (see below for a detailed discussion). In essence, VMSs are only allowed to 
form out of gas that has never perviously experienced star formation, in merging DM halos 
that cross M crit (z) for the first time. 

We only consider mergers in which the two progenitor halos are within a factor of 100 
of each other in mass, although our results are quite similar (within a factor of 2) if this 
ratio is as low as 10 or as high as 300. We impose this limit because very minor mergers will 
not set off a burst of VMS formation, and also because they take a long time. Dynamical 
friction cannot cause halos to merge within a Hubble time if their mass ratio is large. Any 
uncertainty in the allowed mass ratio is significantly smaller than the uncertainty in the 
value of 77*. 

Figure 1 shows the resulting star formation histories. We here presume that VMS 
formation will have ceased at latest by z ~ 7, close to the estimated redshift of reionization 
(Barkana 2002, and references therein). The precise redshift where VMS formation ends is 
uncertain, and is determined by the feedback effects that we will discuss below. For the 
sake of comparison, in Figure 1 we also show an analytical fit to the star formation history 
obtained by Springel & Hernquist (2002) from a series of cosmological simulations that 
incorporate a self-regulated model of quiescent star formation in a multi-phase interstellar 
medium. This latter history describes the formation of a normal, predominantly low-mass 
population of stars (see §3). Note that the Springel & Hernquist results were obtained for a 
cosmology with a slightly smaller value of Qb = 0.04 than for our current analysis. 

2.2. Radiative Feedback and the First Transition in the Star Formation 

History 

In a hierarchical model of structure formation, the very first stars will form in small DM 
halos where gas cooling has to rely on the presence of molecular hydrogen. Initially, therefore, 
star formation proceeds along the molecular branch in Figure 1, and is characterized by a 
top-heavy IMF. This latter prediction follows from the microphysics of H 2 cooling, leading 
to characteristic values of the gas temperature and density, thus imprinting a characteristic 
mass scale of a few 100 M Q (Bromm et al. 2002). The very first stars derive from 3 — 4<r 
peaks and, albeit rare, efficiently produce Lyman- Werner (LW) photons with energies (11.2- 
13.6eV) just below the hydrogen ionization threshold. The LW photons are therefore free 
to propagate through an otherwise neutral IGM, readily destroying H2 in neighboring halos 
(Haiman, Rees & Loeb 1997; Haiman, Abel & Rees 2000; Ricotti, Gnedin & Shull 2001; but 



6 



see Ciardi, Ferrara & Abel 2000a). 

Very massive stars produce approximately ~ 10 4 LW photons per stellar baryon (Bromm 
et al. 2001b). Roughly, one can then estimate the redshift where star formation in low-mass 
halos (relying on H 2 cooling) is no longer possible by demanding that there be of order one 
dissociating LW photon for each molecule in the IGM: 



Here, pe = 4.1 x 10~ 31 g cm' 3 is the comoving density in baryons, /lw the escape fraction 
of LW photons from the star-forming halo, /diss the fraction of absorptions that lead to 
dissociation, / H2 the molecule fraction in the IGM, and 



is the comoving density in stars, with ip*(z) being the comoving star formation rate (SFR), 
in this case given by the molecular track in Figure 1, and the other symbols have their usual 
meaning (e.g., Bromm & Loeb 2002). 

In Figure 2, we evaluate this criterion with /lw — 0.1, /diss — 0.1 (Glover & Brand 

2001) and /h 2 = 10~ 3 , finding z±^2 > 30. We interpret this redshift as the approximate 
boundary between the first two epochs of cosmic star formation. Notice that this is a rather 
conservative estimate, as the fraction of gas in molecular form is only ~ 10~ 6 — 10~ 3 , and 
the required number of LW photons might therefore be significantly smaller than the value 
derived from equation (2). It is also worth pointing out that such an early transition redshift 
is a direct consequence of the predicted top-heavy IMF for the first stars, with a factor > 10 
enhancement in the production of LW photons per stellar baryon compared to the case of a 
standard Salpeter IMF (Bromm et al. 2001b). We have verified that this transition redshift 
is not very sensitive to variations within factors of a few in /lw and /diss- 

At lower redshifts, star formation can now only proceed in halos that are massive enough 
to cool via atomic hydrogen lines (corresponding to masses M > 10 8 M ). To form stars, the 
gas has to cool below the ~ 10 4 K that can be reached by atomic line cooling. Ultimately, 
therefore, star formation in these massive halos still relies on cooling due to H_2 which is 
expected to form in dense regions that are shielded from the LW background (Oh & Haiman 

2002) . The same physical reasoning that leads to the prediction of a top-heavy IMF in small, 
~ 10 6 M , halos should then also apply to these more massive halos, regardless of whether 
H 2 is initially present in them or not. The formation of stars along the atomic track (see 
Figure 2) during epoch 2 consequently proceeds with the same top-heavy IMF as in epoch 
1. 



10 /Lw/dissP*(^l-2) ~ /h 2 PB • 



(2) 




(3) 
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Fig. 2. — The transition from epoch 1 to 2. Comoving star formation rate (SFR) in units 
of M yr~ 1 Mpc" 3 function of redshift. Dotted lines: VMS formation rate calculated for 
cooling due to atomic and to molecular hydrogen, respectively (from Figure 1). Solid line: 
The actually realized star formation history. Once LW-photodissociation feedback becomes 
effective at z ~ 30, the SFR jumps from the molecular to the atomic track. This jump marks 
the transition from the first to the second epoch of high-z star formation. 
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The break in the star formation rate at z = 30 is artificially sharp due to the simplicity 
of our modeling. In reality, this is likely to be a smoother transition because the molecular 
gas will not be dissociated instantaneously, and will not be completely dissociated in more 
massive halos (Machacek, Bryan & Abel 2001). We do believe, however, that this transition 
will be rapid. The mean free path of LW photons at z = 30 is A m f p ~ 10 Mpc (physical), 
much larger than the typical separation of 10 6 M halos at this redshift (~ 30 kpc). Photons 
travel this mean free path in ~ 10 7 yr, which is short compared to the Hubble time. For 
these reasons, a roughly uniform background of LW photons should be set up very quickly. 
The number of sources increases exponentially with time, so the background radiation level 
will increase rapidly until it is intense enough to dissociate most of the molecules in the 
universe. 

The timescale for photodissociation of a molecule is inversely proportional to the inten- 
sity of the radiation field in the Lyman- Werner bands, Jlw, and also depends on how well 
the molecule is shielded from the radiation. Oh & Haiman (2002) give the relationship as 

t diss = 2.1xl0 4 J 2 - 1 V s lL d yr, (4) 

where J21 = Jlw/10 -21 erg cm -2 s _1 Hz -1 sr _1 , and /shield takes into account the effects of 
H 2 self-shielding and is given by /shield = min[l, (iVH 2 /10 14 cm~ 2 )~ - 75 ] (Draine & Bertoldi 
1996). Thus, the dissociation time at fixed column density depends only on the background 
radiation intensity, approximately given by 

J\JN — ] Q ^ C p*(^) ex P[- r IGM] , (5) 

47rm H 

assuming 10 4 LW photons are produced per baryon of star formation (see above). The 
constants in this expression have their usual meaning. The optical depth of the IGM to 
LW photons is estimated to be tigm — 2 — 3 (e.g., Ciardi et al. 2000a; Ricotti et al. 2001), 
reducing the UV flux by about an order of magnitude. Our formula for Jlw neglects the 
redshifting of photons out of the LW bands, but photons travel a mean free path in a small 
fraction of the Hubble time and so should be absorbed long before redshift effects become 
important. 

In Figure 3 we show the dissociation time compared to the Hubble time as a function 
of redshift for different neutral hydrogen column densities. We also show the lifetime of a 
VMS for comparison. We have calculated the curves in Figure 3 for a molecular fraction 
of 10~ 3 , and with r IGM = 2. A 10 6 M Q halo at z = 30 would have a hydrogen column 
density of ~ 4 x 10 20 cm -2 if it were a uniform sphere, and about 10 times larger if it had 
an isothermal density distribution (Glover & Brand 2001). To relate this to the rarity of 
the peaks, we again consider z = 30. A 3-sigma peak has a mass of only 4 x 10 3 M Q , and 
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Fig. 3. — The photodissociation time for molecular hydrogen as a function of redshift. The 
heavy solid line is the Hubble time, the light solid line is the lifetime of a VMS (~ 3Myr), 
and the four broken lines are the dissociation times in halos with different neutral hydrogen 
column densities. Dissociation times decrease towards lower redshifts because the soft-UV 
background radiation is increasing rapidly in intensity. 
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so does not collapse. The extremely rare peaks in the density field are thus responsible 
for dissociating the hydrogen molecules. The masses of halos corresponding to 4-, 5-, and 
6-sigma peaks are ~ 5 x 10 5 , 1 x 10 7 , and 1 x 10 8 M Q , respectively. For a truncated 
isothermal sphere, this corresponds to column densities towards the center of the halos of 
Nr ~ 6 x 10 21 , 2 x 10 22 , and 4 x 10 22 cm -2 . Comparing these column densities with the 
respective curves in Figure 3, we see that the dissociation timescales are always less than 
the Hubble time for even the rarest objects by z = 30. Thus we tentatively conclude that 
this transition is indeed likely to be rapid. 

2.3. Chemical Feedback and the Second Transition in the Star Formation 

History 

Here, we make the simplifying assumption that all VMSs have masses in the range 
140 < M < 260 M Q , and that they consequently explode as pair- instability supernovae 
(PISNe) which completely disrupt the star without leaving behind a compact remnant (Fryer 
et al. 2001; but see Umeda & Nomoto 2002). PISNe are predicted to have substantial metal 
yields, of order y = M z /M* ~ 0.50 (Heger & Woosley 2002). The heavy elements that are 
produced in the SN explosions will enrich the host halo, and a fraction, / m i x , of them will 
escape into the general IGM as a result of SN-driven outflows (e.g., Madau , Ferrara & Rees 
2001; Mori, Ferrara & Madau 2002; Thacker, Scannapieco & Davis 2002). While in §3 we 
predict that some low mass star formation will occur simultaneously with VMS formation, 
these stars will explode as type II supernovae which have a much lower metal yield (~ 1M 
per supernova). This means that their contribution to the metallicity of the universe is small 
as long as there is comparable star formation in VMSs, and we thus ignore any contribution 
from low mass stars in the following calculation. 

With the star formation history described in Figure 2, we now calculate the resulting 
metal enrichment of the IGM 

ZlGM(z) = yimix^-^ , (6) 
PB 

where p*(z) is the cumulative density in stars as given in equation (3). The crucial uncer- 
tainty is the fraction, / m i x , of metals that are able to escape the star forming system, and 
to get mixed into the general IGM. For simplicity, we assume that all the metals escaping a 
star forming halo are mixed uniformly into the IGM. To bracket the likely cases, we calcu- 
late Ziq M (z) for two values of the mixing efficiency: f mix = 0.05 and 0.5. It appears likely 
that the presence of metals is the key determinant in ending the epoch of high-mass star 
formation (Omukai 2000; Bromm et al. 2001a). These authors find that very massive stars 
can no longer form once the gas is enriched to a level in excess of Z ciit ~ 1O~ 3 - 5 Z . The 
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criterion for the second transition in the cosmic star formation history is then 

Zigm( z 2^3) ~ Z crit . (7) 

In Figure 4, we show that this condition results in the approximate range, 15 < £2^3 ^ 20, 
with the limits corresponding to the two assumptions for the mixing efficiency. In a somewhat 
different context, Schneider et al. (2002) have carried out a similar calculation, finding a 
transition redshift, z > 20, that is broadly consistent with our estimated range. 

In reality, not all locations in the universe will have a metallicity close to the mean 
value calculated in equation (6). In fact, the gas in a star forming halo is likely to be rapidly 
enriched to a level ~ Z CT - lt locally (see §3.3), in contrast to the rather slow enrichment of the 
general IGM. A substantial fraction of the Pop III metal production, however, is expected 
to escape from the shallow DM potential wells at high redshifts (e.g., Madau et al. 2001; 
Mori et al. 2002). In addition, the enriched gas has to travel much shorter distances between 
neigboring halos at these early times, and it might therefore have been easier to establish 
a uniform metal distribution in the IGM. The mean IGM metallicity in equation (6) would 
then indeed be a good indicator of when the universe as a whole will undergo the transition 
from predominantly Pop III to predominantly Pop II star formation. Clearly, a more realistic 
understanding of how the metals are mixed, and to what degree of homogeneity, can only 
be gained from numerical simulations. 

In addition to the metal enrichment due to supernovae, the IGM could also be enriched 
due to winds from the first stars. The importance of such winds, however, is predicted to 
be small in the absence of metals and dust grains in the stellar atmosphere (e.g., Kudritzki 
2002). 

2.4. Implications 

The sequence of star-forming modes proposed in this section is expected to apply equally 
to an individual halo (going through different 'stages' of star formation), and in a statistical 
sense to the universe as a whole (going through different 'epochs'). It will be very interesting 
to determine the degree of synchronization in the star-formation modes between spatially 
disjoint regions of the universe. Is it possible that stars in a given location form in a different 
mode from that in a neighboring region? Or is the VMS feedback sudden enough to coor- 
dinate the star formation activity in a global fashion? We emphasize that answers to these 
important questions must await the results from full numerical simulations, and cannot be 
obtained from the simple analytical model presented in this paper. 

By examining Figure 4, it appears that the IGM could reach a level of metal enrichment 
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z 

Fig. 4. — The transition from epoch 2 to 3. IGM metallicity in units of Z as a function of 
redshift. The metallicity is calculated with the star formation history in Fig. 1, assuming 
a heavy element yield of 50%, and two different mixing efficiencies. Solid line: Metal en- 
richment for / m i x = 0.05. Dotted line: Metal enrichment for f mix = 0.5. Horizontal dashed 
line: The critical metallicity above which massive stars are no longer able to form. Vertical 
dot- dashed lines: The redshift at which the IGM metal enrichment reaches Z cr i t , marking 
the transition from the second to the third epoch of star formation at high z. These two 
lines bracket the plausible range, 15 < z < 20, where the transition is likely to occur. 
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from VMS formation which is already close to that inferred for the low-column density 
Lya forest at z ~ 5 (Songaila 2001). Within the context of our model, the IGM would 
have been enriched with a 'bedrock' of metals, corresponding to Z\ GM ~ Z crit , already at 
-22^3 ^ 15. Again, a more realistic assessment of this hypothesis requires sophisticated 
numerical simulations. 



3. Stellar Populations at High z 

3.1. Three Basic Populations 

As a consequence of our picture for high redshift star formation, we predict three distinct 
stellar populations at z > 5 (see Table 1). VMSs form out of gas with Z < Z crit and constitute 
the 'classic' Population III. Low-mass stars, with typical masses of ~ 1 M , that form out 
of already enriched gas with Z > Z crit make up Population II. These stars are actually 
observed, e.g., in the halo of our Galaxy. Finally, we propose a mechanism for producing a 
novel Population II. 5: intermediate-mass stars, forming out of gas with Z < Z CT[t , that owe 
their existence to the mechanical, SN-driven feedback from VMSs. 

3.2. The Formation of Population II. 5 Stars 

We assume that the immediate progenitor of a star is a centrally concentrated, self- 
gravitating gaseous clump with a mass close to the Bonnor-Ebert value (e.g., Palla 2002) 

/ T \ 3/2 / n \ -V2 

M BE ^ 700 M ( — ) . (8) 

° V200K7 V10 4 cm- 3 / v ; 

Here, the temperature and density are normalized to the characteristic values that follow 
from the microphysics of H 2 cooling (see Bromm et al. 2002). Depending on the details of 
how the gas is accreted onto the nascent hydrostatic core in the center of the gravitationally 
unstable clump (e.g., Omukai & Palla 2001; Ripamonti et al. 2002), the resulting stellar mass 
is expected to be somewhat smaller. We take this uncertainty into account by expressing 
the final stellar mass as M* ~ aM BE , and choose the efficiency to be a ~ 0.3. This value is 
close to that inferred for the formation of stars in the present-day universe (McKee & Tan 
2002). For Pop III stars we then have M* > 200 M , in agreement with the results from 
numerical simulations. Notice that the effective star formation efficiency introduced in § 2.1 
is given by i]* ~ f CO oi a , where / coo i is the fraction of the gas that is incorporated into the 
cold, dense clumps. 
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As pointed out in § 2.3, the primordial gas cannot cool below ~ 200 K as long as 
Z < Z ci[t . A decrease in the stellar mass, therefore, has to rely on achieving higher densities, 
beyond the value of ~ 10 4 cm -3 that characterizes Pop III star formation. We propose that 
such a boost in density can naturally occur in the cooled, dense shell behind a SN-driven 
shock wave, in the following way. 

An exploding VMS drives a blast-wave into the surrounding medium that shock-heats 
the gas to a post-shock temperature approximately given by T ps ~ (mH/fe)^, where w s h is 
the (time-dependent) shock speed. Initially, the shock will evolve adiabatically with only a 
modest increase in density, n ps /no < (7 + l)/(7 — 1) =4 for 7 = 5/3. Here, no is the density 
of the undisturbed, preshock gas. Eventually, however, the shocked gas will be able to cool 
and, evolving roughly isobarically (e.g., Shapiro & Kang 1987; Yamada & Nishi 1998), reach 
much higher compression factors. 

To assess whether the post-shock gas can cool sufficiently, we have implemented a simple 
one-zone model, solving for the thermal and chemical evolution in a homogeneous parcel of 
gas. The chemical reaction network and the cooling processes pertaining to primordial gas 
are described in Bromm et al. (2002). Under the assumption of isobaricity, the density evolves 
according to n ~ n T ps /T. Assuming that the post-shock gas is initially fully ionized, the 
only free model parameters are uq and u s h- In Figure 5, we show the resulting temperature 
evolution for uq = 10 2 cm -3 , characteristic of the gas in the vicinity of a VMS (see Bromm 
et al. 2002), and u s ^ = 200 km s _1 , corresponding to T ps ~ 5 x 10 6 K. In this case, the gas 
cools to Tf ~ 200 K within t < 10~ 2 t#(,2 = 20) ~ 10 6 yr, which is shorter than the lifetime 
of a VMS. The overall compression in the case of Figure 5 is then rif ~ n T ps /Tf ~ 10 4 no- 
In general, the characteristic mass of a Pop II. 5 star can be written as 

M,(II.5) - 10 M ( "° \ " /2 ( Ush T ) " . (9) 
* K ; °Vl0 2 cm- 3 7 V200kms-V 

This expression implicitly assumes that the post-shock gas is able to cool to Tf ~ 200 K in 
a time that is short compared to the Hubble time. This requirement constrains the shock 
speed to be u s h < SOOkms -1 . The characteristic mass of a Pop II. 5 star is consequently 
reduced by at least an order of magnitude compared to the Pop III case. 

Within our model, the formation of Pop II. 5 stars is directly linked to the death of a 
Pop III VMS. An important quantity then is the efficiency, 77, of the formation process which 
we define as the mass ratio of intermediate mass (Pop II. 5) stars to very massive (Pop III) 
stars 

M„. 5 eAM shell 

^^~MUmp (10) 

where AM shell is the mass of the cooled, dense shell, e the efficiency of forming stars out of 



this material, and M*(III) ~ 250 M Q the typical mass of a Pop III star. We assume that 
roughly one-half of the mass that has been swept up by the end of the energy-conserving 
(Sedov- Taylor) phase ends up in the dense shell 



AM, 



shell 




(11) 



where E ~ 10 53 erg is the kinetic energy of a PISN (Fryer et al. 2001). Taking the shock 
velocity at the onset of the radiative phase to be u sh ~ 300 km s _1 , we thus estimate 
AMghcii ~ 10 4 ' 5 M©. In equation (10), we conservatively assume that only a small fraction of 
the dense material in the shell is available to form stars, e ~ 0.01, and we therefore predict 
the efficiency of forming Pop II. 5 stars to be small: rj ~ 10 2 e ~ I. 

We would like to stress that the considerations presented here can only serve to physi- 
cally motivate the possibility of such an intermediate stellar population. To fully explore the 
viability of our proposal, we plan to investigate the propagation of a VMS SN with numerical 
simulations (for numerical work on shock-induced star formation in the present-day universe, 
see Vanhala & Cameron 1998 and references therein). 



We propose that the three stellar populations discussed here can form only in certain of 
the epochs introduced in § 2. According to our metallicity criterion, Pop III star formation 
occurs exclusively during the first two epochs. As we have argued above, the formation of 
Pop II. 5 stars is induced by the death of a Pop III star. Only the second epoch, however, is 
conducive to such a triggered star formation mechanism. In epoch 1, Pop III stars form in 
low-mass halos that are already severely disrupted by a single PISN. It is only later, in epoch 
2 of the cosmic star formation history, that the then substantially more massive halos can 
survive a VMS explosion, and consequently enable the triggered formation of intermediate- 
mass, Pop II. 5 stars. To see this more clearly, consider the binding energy of a halo of mass 
M, collapsing at redshift z (Barkana & Loeb 2001): 



Using equation (1), we evaluate this expression for M ~ M crit (z ~ 25), resulting in 



3.3. Implications 
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t/t H 

Fig. 5. — Thermal evolution in shocked gas. Gas temperature vs. time (in units of the 
Hubble time at z ~ 20). Here, the shock has a velocity of u sh = 200 km s" 1 , and propagates 
into a medium with density n = 10 2 cm~ 3 . The initial, very sudden decrease in temperature 
is due to the onset of atomic (He and H) line cooling. The subsequent drop is brought about 
by molecular (H2) cooling. It is evident that the shocked gas is able to cool to ~ 200 K 
within a small fraction of the Hubble time at z — 20. 
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It is evident that the criterion for the formation of Pop II. 5 stars, E b > E ~ 10 53 erg, is 
fulfilled in halos that cool via atomic hydrogen lines, but not in those that cool via H 2 . 

Within the framework presented in this paper, Pop II stars can clearly form in epoch 
3, at which point we assume that most of the IGM is enriched to a level in excess of Z crit . 
It is likely, however, that Pop II star formation already ensues during epoch 2. The metal 
enrichment due to VMSs is so efficient that the host system is expected to reach the critical 
metallicity well before the general IGM does. In fact, one can estimate that already one 
single VMS could enrich the ~ 10 M of gas in the typical star-forming halo in epoch 2 to 
the critical level: 

7 M z 

^halo 



M 



B 



10 2 M, 



10 3 ' 3 Z > Z cr it , 



- 10 7 M 

where we have assumed for simplicity that none of the metals have escaped the halo, corre- 
sponding to the limiting case of / m i x = (see § 2.3). 

During epoch 2, therefore, we expect all three stellar populations to occur almost simul- 
taneously in a given star forming system. In Table 1, we summarize the association of the 
various stellar populations with the distinct epochs of star formation. 

We now tentatively suggest how these stellar populations might relate to the different 
classes of SNe as proposed by Qian & Wasserburg (2002). These authors (henceforth QW) 
account for the observed abundance patterns in low-metallicity Galactic halo stars within 
the framework of three distinct classes of SNe. At the beginning of the nucleosynthetic chain, 
a generation of very massive stars explodes to enrich the primordial material mostly with Fe, 
but producing no elements beyond the iron peak. It seems natural to identify the QW VMS 
SNe with the Pop III PISNe discussed here. According to QW, at some later stage in the 
chemical enrichment history, more normal, type II, SNe begin to contribute nucleosynthetic 
products. More specifically, QW distinguish between a class of high-frequency and low- 
frequency type II SNe, named SN II(-f/) and SN ll(L), respectively. The H events produce 
mainly the heavy r-process elements but no Fe, whereas the L events contribute mainly 
the light r-process elements as well as Fe. These two classes of SN II are postulated on 
purely phenomenological grounds, and it is presently not known what kind of progenitor 
stars could physically give rise to them. There are, however, some indications that the H 
events correspond to the explosion of comparatively low-mass stars that are only barely able 
to reach the SN stage (see QW and references therein). Given that Pop II. 5 stars are more 
massive than Pop II stars by ~one order of magnitude, we speculate that Population II. 5 
might predominantly give rise to SNe H(L), whereas Population II could lead to both the H 
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and L events. Again, we summarize these proposed identifications in Table 1. 

Regardless of the detailed correspondence discussed above, our general model for the 
character of star formation at high redshift seems to provide a promising astrophysical con- 
text for the QW proposal. More work, both observational and theoretical, is required to 
explore this possible connection further. 



4. Observational Consequences 
4.1. Supernova Rates as Star Formation Tracer 

4-1.1. Modeling 

The cosmic supernova rate as a function of redshift is potentially an excellent tracer of 
the different modes of star formation envisaged in this paper. A number of authors have 
previously advocated this method of tracing star formation, predicting supernova rates based 
on the metallicity of the IGM at high redshifts (Miralda-Escude & Rees 1997), or on the 
empirically determined star formation rate out to z ~ 5 (e.g., Madau, Delia Valle & Panagia 
1998; Sadat et al. 1998; Dahlen & Frannson 1999). In our study, we use theoretical estimates 
of the cosmic star formation rate to predict the number of supernovae from VMSs and from 
regular stars (Pop II) out to much higher redshifts. Observational determinations of these 
rates will constrain the three epochs of star formation which we are proposing. Including 
observational selection functions is beyond the scope of this paper. We will, however, present 
rough estimates of the flux limit and wavelength coverage needed in order to observe SN 
explosions at high redshift. The observational capability of NGST has recently been studied 
for high- 2: type II SNe (Dahlen & Frannson 1999), and in a preliminary fashion for PISNe 
from VMSs (Heger et al. 2002). 

We use the results of recent cosmological simulations by Springel & Hernquist (2002) 
to model the Pop II star formation history of the universe. These simulations agree roughly 
with observations of the star formation rate out to z ~ 5, but extend to much higher redshifts 
than are currently probed by direct observations, and also reproduce the mean density of 
stars in the universe. These authors provide a simple analytical fit to the SFR as follows: 

^{z',U) = K- ■ =- -, 13 

b — a + aexp[b(z — z m )\ 

where the normalization is K = 0.15 M yr _1 Mpc~ 3 and the fitting parameters are a = 3/5, 
b = 14/15, and z m = 5.4. The time delay between the formation and explosion of a massive 
star is < 20 Myr, and therefore the Type II supernova rate should trace the star formation 
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rate almost exactly. For a Salpeter IMF, there is one Type II supernova (SN II) per ~ 150 M 
of star formation, and we consequently divide the SFR by this factor (which we call Msn(H) ) 
to obtain the Pop II SN rate. Hernquist & Springel (2002) find that the SFR at z < z m in 
their simulations actually has a power-law dependence on the expansion rate, rather than 
evolving exponentially with redshift; however the fit given in the above equation is sufficient 
for our present purposes. 

For the PISN rate we use the VMS star formation rate in Figure 2. The typical lifetime 
of a VMS is ~ 3 Myr so we can again assume the SN rate traces the SFR closely. We take 
typical VMSs to have masses of 250 M and assume they all die as supernovae to calculate 
the VMS supernova rate (i.e. Msn(HI) = 250 M ). VMSs are predicted to end their lives as 
pair-instability supernovae if M < 260 M , and these should be very energetic and luminous 
events. Bear in mind, however, that the mass function, and hence end state of Pop III stars 
is still somewhat uncertain, making this supernova rate a crude estimate. 

We also calculate the rate of Pop II. 5 supernovae, using the result from the last section 
that the Pop II. 5 SFR is directly proportional to the Pop III SFR, with constant of pro- 
portionality 7] ~ 1. Due to the uncertain nature of the Pop II. 5 IMF we must make simple 
assumptions to get the mass in stars formed per supernova. The characteristic stellar mass 
was shown to be about 10 M , so we assume that half of the stellar mass goes into stars that 
can explode. We further assume that the average mass of a Pop II. 5 supernova progenitor 
is ~ 20 M . From this we find that Msn(H-5) ~ 40 M for Pop II. 5 supernovae. We expect 
these supernovae to be of Type II, although with very low metallicity. 

We convert these rates per unit volume into rates observed per unit time and per unit 
redshift interval: 



Here dV is a comoving volume element, Msn is the mass in stars formed per supernova for 
the relevant stellar population, f sky the fraction of the sky observed, and c/H the Hubble 
length. Note that the Pop II. 5 and Pop III supernova rates differ only in their respective 
values for M* calculated above. The respective time intervals in the observer and emitter 
frames, dt Q b s and dt em , are related by the factor of (1 + z) in the denominator. The comoving 
distance to redshift z is 



dN 



dN dV 

(1 + z)dt em dV dz 
1 ip^z) c 47r f s k y r 2 {z) 



dt bsdz 



M SN (1 + z) H ^Q m (i + z f + Q A ' 



(14) 
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4-1.2. Results 

The results for our cosmological parameters are shown in Figure 6, where we plot the 
number of supernovae that could be observed per year and per square degree. Note that this 
figure does not incorporate magnitude limits or observational selection functions, but does 
simply show the predicted number of supernovae that can in principle be observed. 

The redshift distribution of SNe follows the star formation rate, peaking at z ~ 5 for 
Pop II, and showing the same break at z = 30 for PISNe as in Figure 2. We expect the 
transition from predominantly PISNe to Pop II SNe to occur around a redshift of z ~ 15 — 20 
(see Figure 4), and we therefore show the rates for both types of supernova in this redshift 
range. As this transition occurs, we expect the Pop II SN rate to gradually drop below the 
solid line towards higher redshifts as Pop II star formation switches off. Similarly we expect 
the Pop III (Pop II. 5) rate to drop below the dashed (dotted) line towards lower redshifts 
as Pop III star formation ceases. We do not expect any VMSs at z < 15, but we extend our 
calculation down to z = 7 because the reionization of the universe at roughly this redshift 
provides a robust lower limit for the end of VMS formation (e.g., Oh et al. 2001). That 
the PISN and Pop II SN rates are so similar at the transition redshifts is because the star 
formation rate for Pop III at z ~ 15 which we predict analytically is similar to the rate (for 
Pop II) found in the simulations of Springel & Hernquist (see Figure 1). Again, we emphasize 
that the break in the VMS supernova rate at z = 30 is artificially sharp, reflecting the abrupt 
transition in the underlying star formation rate. 

In terms of total numbers of supernovae, we predict ~ 2 x 10 6 PISNe per year over 
the whole sky (above z = 15), or about 50 per square degree per year. This still comprises 
only ~ 0.4% of all supernovae, as we find ~ 4 x 10 8 SNell per year (all sky), a value 
comparable to previous estimates (Miralda-Escude & Rees 1997; Madau et al. 1998). This 
number can be understood by a simple order of magnitude estimate as follows. About 10% 
of the baryons in the universe are in stars by the present time, giving a local density in stars 
of p* ~ 0.1f2bp cr it ~ 6 x 10 8 M /Mpc 3 for our cosmological parameters (Peebles 1993). The 
size of the observable universe is ~ 10 Gpc out to z = 5 (most of the star formation happens 
later than z = 5 because there is very little physical time before this), so the total mass in 
stars in this volume, ignoring redshift evolution effects, is ~ 10 21 - 5 M Q . About 10% of the 
mass in stars goes into stars which can explode, so the mass in supernovae is ~ 10 20 5 M . 
The mass weighted average for a supernova progenitor's mass is ~ 30 M Q , and the age of 
the universe is 14 Gyr. The number of supernovae per year is then 

d/V 10 20 - 5 

-P- = 7 x 10 8 SNe yr" 1 , (16) 

dt obs 30 x 14 x 10 9 yr 3 K ' 

similar to our numerical integration of the calculated SN rate. All calculations of the SN 
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Fig. 6. — Supernova rate as a function of redshift for Type II and VMS supernovae. Pop II 
and Pop II. 5 refer to Type II supernovae from Pop II and Pop II. 5 stars, respectively. Pop III 
refers to VMS supernovae. The assumed star formation rates are detailed in the text. We 
show the number of supernovae per observed year in one square degree of the sky, plotted 
as a function of redshift. 
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rate based on the cosmic SFR are expected to yield an answer around this value. 



4-1.3. Observational Prospects 

We now address the question of whether any of these supernovae at very high redshift 
could be observed in the foreseeable future. Type la supernovae have a peak B-magnitude 
of M B ~ —19, with Type II supernovae being about 2.5 magnitudes fainter (e.g., Filippenko 
1997). VMSs are predicted to end their lives as pair instability supernovae if their mass is 
in the range 140 M < M < 260 M (Heger & Woosley 2002), but they have not yet been 
observed and their luminosity is not well known. Heger & Woosley (2002) have shown that 
they can produce large masses of 56 Ni during supernova nucleosynthesis. If the light curve is 
powered by 56 Ni decay, then these supernovae should be brighter than Type la's by a large 
factor (the ratio of their 56 Ni production). Preliminary results from Heger et al. (2002) show 
that reality is somewhat more complicated, and these authors find that a PISN is only a 
factor of a few brighter than a type la. Because of this uncertainty in luminosity, we will 
calculate the apparent brightness of a type la SN at large redshifts, and scale this up for the 
brighter PISN. The apparent brightness in a flat universe is given by 

'•<*> = MTT$W> • (17) 

where L v is the monochromatic luminosity of the source at redshift z. For z = 20, Mb = —19, 
and our cosmological parameters we obtain 

/„(A obs ~ 9/im, z = 20) ~ 5 nJy . (18) 

The observed wavelength of 9/xm is the redshifted 5-band. If the flux does not drop off 
too much in the rest frame near-UV, then there should be a similar flux level in the near- 
IR (A < 5/im) where NGST should have sub-nJy sensitivity (Barkana & Loeb 2001). This 
indicates that even a Type la SN should be observable with NGST out to z ~ 20, and PISNe 
should be easily observable if they are at least a few times brighter. Type II supernovae, on 
the other hand, will be very hard to detect at these high redshifts because they are typically 
a factor of 10 fainter. For this reason, even though we predict Pop II. 5 supernovae to be 
more numerous than PISNe, we expect that they will not be nearly as useful as observational 
probes of the high redshift universe. Of course observing them would be useful for assessing 
the efficiency of low mass star formation, and for testing our shock compression model. 

Marri & Ferrara (1998) have discussed the feasibility of using gravitational lensing mag- 
nification due to intervening mass concentrations to observe high redshift supernovae more 
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easily, and potentially bring even type II supernovae within reach of NGST at z > 10. Sim- 
ilar considerations also apply to our calculations, and this may help in obtaining detailed 
observations of PISNe. It is important to be able to identify them reliably, distinguish them 
from other types of supernova, and measure their redshift if we are to use them for cosmo- 
logical purposes. To be able to do this, we need a sample of SNe that are bright enough for 
spectroscopic observations. Since none have been observed locally, we need the additional 
lensing magnification to investigate SNe in detail at high redshift. 



If our mechanism for producing extremely low metallicity, low mass stars (Pop II. 5 in 
our terminology) via shock compression occurs in reality, then some of these stars should 
still be present in the halo of our galaxy. They would of course be very old (~ 14 Gyr), so 
only the lowest mass stars with M < 0.8 M Q would remain (Girardi et al. 2000; Siess, Livio 
& Lattanzio 2002). We can use our model to predict how many such stars we would expect 
to find in a typical Milky Way-sized halo and, given a model for the stellar halo density 
profile, their number density in the solar neighborhood (we use a similar approach to that 
of Hernandez & Ferrara 2001). Significant uncertainties in this calculation are the efficiency 
of producing Pop II. 5 stars (as discussed in §3.2), the fraction of them that end up in the 
halo as opposed to the bulge, and their mass function. 

The star formation rate of Pop II. 5 stars is given by ip*(z; II. 5) = 77^* (2:; III) in the 
redshift range when they form. There is a slight difference to the way ip* (z; III) is calculated 
for this problem, in that we are interested in the stars that end up in the Milky Way (MW) 
halo. Therefore, instead of the usual Press-Schechter mass function in the Santos et al. (2002) 
model, we use the Extended Press-Schechter mass function of halos at redshift z which will 
end up in a halo the mass of the MW at z = (Lacey & Cole 1993). 

For a MW mass (total) of M MW = 10 12 M we find that 10 7 - 5 M of Pop III stars are 
formed in MW progenitor halos, over the redshift range from z = 30 to z = 15 in which we 
expect Pop II. 5 stars to form. Interestingly, because the metal yield of PISNe is ~ 0.5, this 
results in a MW metal enrichment of 
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similar to the metallicity of the lowest metallicity globular clusters and thick disk stars 
([Fe/H]= -2.2, Freeman & Bland-Hawthorn 2002). 



For a given value of 77, this mass in Pop III stars gives us the corresponding mass in 
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Pop II. 5 stars: 

Mll5 , 10 -M e (^)(^). (20) 

We now wish to convert this into the number density of Pop II. 5 stars that could be observable 
in the MW halo today. First, we have argued that the characteristic mass of Pop II. 5 stars is 
~ 10 M , so most of these stars will have died long before the present. We parameterize this 
in a very simple way by assuming that some fraction, j\ m f ~ 0.1 of the mass in stars formed 
goes into stars with M < 0.8 M Q . A second poorly constrained process is where in the MW 
these stars reside. It is expected that most of the earliest forming stars will be incorporated 
into the bulge of the MW, with a small fraction in the halo (White & Springel 2000). We 
introduce /halo ~ 0.1 as the fraction of Pop II. 5 stars which end up in the halo. Finally we 
make the reasonable assumption that most of the Pop II. 5 stars remaining today are near 
the end of their lives (i.e. their mass function is increasing to higher masses at M < 1 M ), 
so we get the number of halo stars by dividing their total mass by 0.8 M Q . This gives us the 
total number of Pop II. 5 stars in the MW halo: 

Nn. 5 ~ ID 5 5 M (X) . (21) 

To relate this to the observed number density of low metallicity halo stars we distribute the 
Pop II. 5 stars in a halo with density profile p(r) = p(r )r /(rp +r 3 ), where the solar orbital 
radius, r = 8.0 kpc (Reid 1993), is assumed much larger than the core radius, r c , of the 
halo profile (needed for finite halo mass). This distribution only diverges logarithmically at 
large radii, so the mass of the halo is quite insensitive to the uncertain cutoff radius and core 
radius (we take r cut /r core = 300 for deflniteness) . Using this profile, together with the total 
number of stars from the previous paragraph, the number density of Pop II. 5 stars in the 
solar neighborhood is 

— 10st *(^)(Tj)(tf)- < 22 > 

We can make a crude comparison of this prediction with the results of observational 
searches for metal poor halo stars (e.g., Beers et al. 2000). Using their Figure 5, we see 
that most of their stars are within 2 kpc of the sun, which we take to be the radius of their 
survey (covering 1940 square degrees of the sky). Also from that figure, they have ~ 5 stars 
with [Fe/H]< —3.5, and none with [Fe/H]< —4.0. If we identify these stars with Pop II. 5, 
this gives them an observed space density of n ~ 3 kpc -3 . This is similar to our theoretical 
estimate above, given the uncertainty associated with some of our parameters. 

The observed stars all have a metallicity very close to our critical metallicity, so it is 
not clear if they are very low metallicity Pop II stars or if they were formed by our Pop II. 5 
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mechanism. We have shown that, with reasonable modeling assumptions, it is possible to 
produce these stars from primordial gas, and that their metallicity could be due to pollution 
from Pop III supernovae which triggered their formation (it could also be partly due to 
pollution from swept up gas as they travel through the galactic disk) . Current observational 
searches are beginning to give us a clearer picture of the number of very low metallicity stars 
in the galactic halo, and we anticipate that this will continue with future observations. On 
the theoretical side, numerical simulations can in principle reduce the uncertainties in our 
Pop II. 5 model by determining the values of 77, /halo, and possibly also of the IMF from the 
mass spectrum of collapsing clumps, although this is a very challenging problem to simulate 
accurately. 



5. Summary and Conclusions 

We have investigated the history of star formation in the high redshift universe, focusing 
on the role played by very massive Pop III stars. We have argued that this history is shaped 
by the various feedback effects exerted by those stars, resulting in three distinct epochs of 
star formation. 

The first impact Pop III stars have on their surroundings is radiative, as their soft 
UV flux dissociates the molecular hydrogen in other minihalos (Haiman et al. 1997, 2000). 
Because of their copious UV emission, and the small molecular fraction (~ 10~ 4 — 10~ 6 ) 
in minihalos, this photodissociation is likely to happen quickly, and fairly completely. A 
numerical investigation of this process was performed by Machacek et al. (2001), who sim- 
ulated this suppression of star formation in low mass halos by incorporating a uniform soft 
UV background into a cosmological simulation. While this assumption does not allow them 
to follow the detailed build-up of the UV radiation field as Pop III star formation switches 
on, it does give a good indication of the overall effect. It was found that the UV radiation 
can effectively suppress star formation due to molecular cooling in low mass halos, delaying 
their collapse from redshift z ~ 30 to z ~ 20. 

It is encouraging that both analytic and numerical calculations agree quite well. In 
Figure 2, we find that once a significant UV radiation field has been set up by z ~ 30, 
further star formation is strongly inhibited until more massive halos start to form, and the 
SFR does not recover to its z ~ 30 value until z < 20. This is all in good agreement with 
the simulation results, except that Machacek et al. (2001) did find that more massive halos 
of M > 10 7 M were able to retain some molecular hydrogen. Hence, these halos were still 
able to cool via molecular lines, indicating that the SFR should not go all the way down 
to the atomic cooling curve, but it is likely to drop most of the way. Future simulations 
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will be able to follow this dissociation process in detail, taking into account the effect of 
self-shielding, with a time-varying radiation field generated by point sources at the sites of 
Pop III star formation. In summary, the radiative feedback from Pop III stars results in the 
transition, at z ~ 30, between the first two epochs of star formation at high redshift. 

The second impact Pop III stars have on surrounding material is chemical, and we 
identify this feedback as the main effect governing the transition from forming predominantly 
very massive stars to forming normal, lower mass ones (e.g., Omukai 2000; Nishi & Tashiro 
2000; Bromm et al. 2001a; Schneider et al. 2002). With the assumption that VMSs can only 
form in gas with a metallicity lower than a critical value, Z CT - lt ~ 10~ 3 ' 5 Z Q , we estimate that 
most of the star-forming material in the IGM has reached this level in the redshift range 
z ~ 15 — 20, which delineates the second transition in the star formation history, between 
epochs 2 and 3. 

An important issue is whether such a well-defined transition redshift between massive 
and more normal star formation actually occurs. This depends on how synchronized different 
parts of the universe are in crossing the critical metallicity threshold, which in turn depends 
on how well the supernova ejecta are mixed through the IGM at high redshifts. It has been 
suggested (Madau et al. 2001) that enrichment becomes more synchronized and uniform at 
higher redshifts. Dark matter halos are of much lower mass and have shallower potential 
wells, making it easier for metals to escape the halo. Supernovae from VMSs are predicted to 
be very energetic, with explosion energies significantly larger than the binding energy of halos 
at z > 20 (Barkana & Loeb 2001). Furthermore, the universe is much denser at high redshift 
and halos are consequently much closer together, again making it easier to pollute the IGM 
uniformly with metals. It is also possible, however, that enrichment at high redshifts is similar 
to that at low redshifts, where it has been shown to be very inhomogeneous and incomplete 
(e.g., Gnedin 1998; Cen & Ostriker 1999; Aguirre et al. 2001a,b). If this were the case, the 
transition redshift at which the mean mass-weighted metallicity of the universe crosses the 
critical threshold would not be a very meaningful quantity. Which of these two cases is true 
is an important and interesting question, and can only be addressed theoretically by more 
detailed cosmological simulations which can resolve star-forming halos at high redshifts. 

The demise of Pop III stars in SN explosions opens up the possibility for a new stellar 
population which we have termed Population II. 5. These stars are characterized by masses, 
M* ~ 10 M Q , intermediate to that of Pop III and Pop II stars, and by extremely low 
metallicities, below the critical value that would normally lead to the formation of VMSs. 
This possibility relies on the gas in the vicinity of the SN being compressed in a radiative 
shock, leading to a significant increase in density, and thus possibly to a reduction of the 
fragment mass, Mre oc p" 1 / 2 , to values much lower than the one realized in primordial 
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gas. During an intermediate epoch, therefore, low- and high-mass stars might form almost 
simultaneously. A similar idea had already been proposed in a prescient paper by Cayrel 
(1986), primarily to explain the so-called G-dwarf problem, i.e., the observed lack of metal- 
poor stars in our Galaxy (see also Kashlinsky & Rees 1983; Nakamura & Umemura 2001). 

It is to be expected that the mass range of Pop II. 5 stars extends down to < 1 M , and 
these stars at the low-mass end of the Pop II. 5 IMF should still be around today. Intriguingly 
and possibly providing a case in point, the recent discovery of the extremely metal poor star 
CS 29498-043 with [Fe/H]= —3.7 and a significant overabundance of Mg and Si, has been 
tentatively interpreted as hinting at a new class of stars (Aoki et al. 2002). Ongoing surveys 
of extremely metal-poor halo stars, much improved in size and quality, should soon be able 
to test our prediction further. 

Many of the uncertainties in our argument have to be addressed with detailed numerical 
simulations, and we plan to do so in future work. The framework presented in this paper 
does provide a coherent context for these numerical studies, and it highlights the important 
physical questions that define the challenge of elucidating the end of the cosmic "dark ages" . 

We are grateful to Andrea Ferrara, Alex Heger, Marc Kamionkowski, Robert Kirshner, 
Avi Loeb, Thomas Matheson, Mike Santos, Volker Springel, and Jerry Wasserburg for helpful 
discussions. We thank the anonymous referee for comments that improved the presentation 
of this paper. This work was supported in part by NSF grants ACI 96-19019, AST 98-02568, 
AST 99-00877, and AST 00-71019. 
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Table 1. Stellar Populations at High Redshifts 





M*/M 


Z/Z Q 


QW SNe a 


Epoch 


Pop III 


> 100 


< -^crit 


VMS 


1,2 


Pop II.5 


• >io 


< -^crit 


SN 11(1) 


2 


Pop II 


>1 


> -^crit 


SN ll(H),(L) 


2,3 



a Tentative identification with the different kinds of SNe proposed by 
Qian & Wasserburg (2002). 



